PAMAM Dendrimer-Drug Interactions: Effect of pH on the Binding and Release Pattern

Understanding the dendrimer-drug interaction is of great importance to design and optimize the dendrimer-based drug delivery system. Using atomistic molecular dynamics (MD) simulations, we have analyzed the release pattern of four ligands (two soluble drugs, namely, salicylic acid (Sal), L-alanine (Ala), and two insoluble drugs, namely, phenylbutazone (Pbz) and primidone (Prim)), which were initially encapsulated inside the ethylenediamine (EDA) cored polyamidoamine (PAMAM) dendrimer using the docking method. We have computed the potential of mean force (PMF) variation with generation 5 (G5)-PAMAM dendrimer complexed with drug molecules using umbrella sampling. From our calculated PMF values, we observe that soluble drugs (Sal and Ala) have lower energy barriers than insoluble drugs (Pbz and Prim). The order of ease of release pattern for these drugs from G5 protonated PAMAM dendrimer was found to be Ala > Sal > Prim > Pbz. In the case of insoluble drugs (Prim and Pbz), because of larger size, we observe much nonpolar contribution, and thus, their larger energy barriers can be reasoned to van der Waals contribution. From the hydrogen bonding analysis of the four PAMAM drug complexes under study, we found intermolecular hydrogen bonding to show less significant contribution to the free energy barrier. Another interesting feature appears while calculating the PMF profile of G5NP (nonprotonated)-PAMAM Pbz and G5NP (nonprotonated)-PAMAM-Sal complex. The PMF was found to be less when the drug is bound to nonprotonated dendrimer compared to the protonated dendrimer. Our results suggest that encapsulation of the drug molecule into the host PAMAM dendrimer should be carried out at higher pH values (near pH 10). When such complex enters the human body, the pH is around 7.4 and at that physiological pH, the dendrimer holds the drug tightly. Hence the release of drug can occur at a controlled rate into the bloodstream. Thus, our findings provide a microscopic picture of the encapsulation and controlled release of drugs in the case of dendrimer-based host-guest systems.

Molecular dynamics simulations of PPI dendrimer-drug complexes

Dendrimeric nanoparticles are potential drug delivery devices which can enhance the solubility of hydrophobic drugs, thus increasing their bioavailability and sustained release action. A quantitative understanding of the dendrimer-drug interactions can give valuable insight into the solubility and release profile of hydrophobic drug molecules in various solvent conditions. Fully atomistic molecular dynamics (MD) simulations have been performed to study the interactions of G5 PPIEDA (G5 ethylenediamine cored poly(propylene imine)) dendrimer and two well known drugs (Famotidine and Indomethacin) at different pH conditions. The study suggested that at low pH the dendrimer-drug complexes are thermodynamically unstable as compared to neutral and high pH conditions. Calculated Potential of Mean Force (PMF) by umbrella sampling showed that the release of drugs from the dendrimer at low pH is spontaneous, median release at neutral pH and slow release at high pH. In addition, Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) binding free energy calculations were also performed at each umbrella sampling window to identify the various energy contributions. To understand the effect of dendrimer chemistry and topology on the solubility and release profile of drugs, this study is extended to explore the solubility and release profile of phenylbutazone drug complexed with G3 poly(amidoamine) and G4 diaminobutane cored PPI dendrimers. The results indicate that the pH-induced conformational changes in dendrimer, ionization states, dendrimer type and pK(a) of the guest molecules influence the free energy barrier and stability of complexation, and thus regulate drug loading, solubility and release.

Intermolecular interactions, charge-density distribution and the electrostatic properties of pyrazinamide anti-TB drug molecule: an experimental and theoretical charge-density study

An experimental charge-density analysis of pyrazinamide (a first line antitubercular drug) was performed using high-resolution X-ray diffraction data (sin theta/lambda)(max) = 1.1 angstrom(-1)] measured at 100 (2) K. The structure was solved by direct methods using SHELXS97 and refined by SHELXL97. The total electron density of the pyrazinamide molecule was modeled using the Hansen-Coppens multipole formalism implemented in the XD software. The topological properties of electron density determined from the experiment were compared with the theoretical results obtained from CRYSTAL09 at the B3LYP/6-31G** level of theory. The crystal structure was stabilized by N-H center dot center dot center dot N and N-H center dot center dot center dot O hydrogen bonds, in which the N3-H3B center dot center dot center dot N1 and N3-H3A center dot center dot center dot O1 interactions form two types of dimers in the crystal. Hirshfeld surface analysis was carried out to analyze the intermolecular interactions. The fingerprint plot reveals that the N center dot center dot center dot H and O center dot center dot center dot H hydrogen-bonding interactions contribute 26.1 and 18.4%, respectively, of the total Hirshfeld surface. The lattice energy of the molecule was calculated using density functional theory (B3LYP) methods with the 6-31G** basis set. The molecular electrostatic potential of the pyrazinamide molecule exhibits extended electronegative regions around O1, N1 and N2. The existence of a negative electrostatic potential (ESP) region just above the upper and lower surfaces of the pyrazine ring confirm the pi-electron cloud.

Cocrystals of Hydrochlorothiazide: Solubility and Diffusion/Permeability Enhancements through Drug-Coformer Interactions

Hydrochlorothiazide (HCT) is a diuretic and a BCS class IV drug with low solubility and low permeability, exhibiting poor oral absorption. The present study attempts to improve the physicochemical properties of the drug using a crystal engineering approach with cocrystals. Such multicomponent crystals of HCT with nicotinic acid (NIC), nicotinamide (NCT), 4-aminobenzoic acid (PABA), succinamide (SAM), and resorcinol (RES) were prepared using liquid-assisted grinding, and their solubilities in pH 7.4 buffer were evaluated. Diffusion and membrane permeability were studied using a Franz diffusion cell. Except for the SAM and NIC cocrystals, all other binary systems exhibited improved solubility. All of the cocrystals showed improved diffusion/membrane permeability compared to that of HCT with the exception of the SAM cocrystal. When the solubility was high, as in the case of PABA, NCT, and RES cocrystals, the flux/permeability dropped slightly. This is in agreement with the expected interplay between solubility and permeability. Improved solubility/permeability is attributed to new drug-coformer interactions. Cocrystals of SAM, however, showed poor solubility and flux This cocrystal contains a primary sulfonamide dimer synthon similar to that of HCT polymorphs, which may be a reason for its unusual behavior. Hirshfeld surface analysis was carried out in all cases to determine whether a correlation exists between cocrystal permeability and drug-coformer interactions.

Structure and interactions of aspirin-like drugs

A pre-requisite for the elucidation of the mechanism of action of aspirin-like drugs, which are believed to exert their pharmacological effects through the inhibition of prostaglandin biosynthesis, is an understanding of their molecular geometry, the non-covalent interactions they are likely to be involved in, and the geometrical and the electronic consequences of such interactions. This has been sought to be achieved through the x-ray analysis of these drug molecules and their crystalline complexes with other suitable molecules. The results obtained from such studies have been discussed in terms of specific typical examples. For instance, antipyrine can form metal and hydrogen-bonded complexes; phenylbutazone can form ionic complexes with basic molecules. Complex formation is accompanied by characteristic changes in the molecular geometry and the electronic structure in both the cases. The results obtained so far appear to indicate that the important common invariant structural features of the fenamates, deduced from crystal structures, are retained even when complexation takes place.

Structural Basis of Drug Resistance by Genetic Variants of HIV Type 1 Clade C Protease from India

Using computer modeling of three-dimensional structures and structural information available on the crystal structures of HIV-1 protease, we investigated the structural effects of mutations, in treatment-naive and treatment-exposed individuals from India and postulated mechanisms of resistance in clade C variants. A large number of models (14) have been generated by computational mutation of the available crystal structures of drug bound proteases. Localized energy minimization was carried out in and around the sites of mutation in order to optimize the geometry of interactions present. Most of the mutations result in structural differences at the flap that favors the semiopen state of the enzyme. Some of the mutations were also found to confer resistance by affecting the geometry of the active site. The E35D mutation affects the flap structure in clade B strains and E35N and E35K mutation, seen in our modeled strains, have a more profound effect. Common polymorphisms at positions 36 and 63 in clade C also affected flap structure. Apart from a few other residues Gln-58, Asn-83, Asn-88, and Gln-92 and their interactions are important for the transition from the closed to the open state. Development of protease inhibitors by structure-based design requires investigation of mechanisms operative for clade C to improve the efficacy of therapy.

Computer-Aided Drug Design of Pyranopyrazoles and Related Compounds for Checkpoint Kinase-1

Checkpoint-1 kinase plays an important role in the G(2)M cell cycle control, therefore its inhibition by small molecules is of great therapeutic interest in oncology. In this paper, we have reported the virtual screening of an in-house library of 2499 pyranopyrazole derivatives against the ATP-binding site of Chk1 kinase using Glide 5.0 program, which resulted in six hits. All these ligands were docked into the site forming most crucial interactions with Cys87, Glu91 and Leu15 residues. From the observed results these ligands are suggested to be potent inhibitors of Chk1 kinase with sufficient scope for further elaboration.

Three-dimensional structure of heat shock protein 90 from Plasmodium falciparum: molecular modelling approach to rational drug design against malaria

We have recently implicated heat shock protein 90 from Plasmodium falciparum (PfHsp90) as a potential drug target against malaria. Using inhibitors specific to the nucleotide binding domain of Hsp90, we have shown potent growth inhibitory effects on development of malarial parasite in human erythrocytes. To gain better understanding of the vital role played by PfHsp90 in parasite growth, we have modeled its three dimensional structure using recently described full length structure of yeast Hsp90. Sequence similarity found between PfHsp90 and yeast Hsp90 allowed us to model the core structure with high confidence. The superimposition of the predicted structure with that of the template yeast Hsp90 structure reveals an RMSD of 3.31 angstrom. The N-terminal and middle domains showed the least RMSD (1.76 angstrom) while the more divergent C-terminus showed a greater RMSD (2.84 angstrom) with respect to the template. The structure shows overall conservation of domains involved in nucleotide binding, ATPase activity, co-chaperone binding as well as inter-subunit interactions. Important co-chaperones known to modulate Hsp90 function in other eukaryotes are conserved in malarial parasite as well. An acidic stretch of amino acids found in the linker region, which is uniquely extended in PfHsp90 could not be modeled in this structure suggesting a flexible conformation. Our results provide a basis to compare the overall structure and functional pathways dependent on PfHsp90 in malarial parasite. Further analysis of differences found between human and parasite Hsp90 may make it possible to design inhibitors targeted specifically against malaria.

Interactions of linear dicationic molecules with lipid A: structural requisites for optimal binding affinity

The structural determinants of the binding affinity of linear dicationic molecules toward lipid A have been examined with respect to the distance between the terminal cationic functions, the basicity, and the type of cationic moieties using a series of spermidine derivatives and pentamidine analogs by fluorescence spectroscopic methods, The presence of two terminal cationic groups corresponds to enhanced affinity, A distinct sigmoidal relationship between the intercationic distance and affinity was observed with a sharp increase at 11 Angstrom, levelling off at about 13 Angstrom. The basicity (pK) and nature of the cationic functions are poor correlates of binding potency, since molecules bearing primary amino, imidazolino, or guanido termini are equipotent, The interaction of pentamidine, a bisamidine drug, with lipid A, characterized in considerable detail employing the putative intermolecular excimerization of the drug, suggests a stoichiometry of 1:1 in the resultant complex, The binding is driven almost exclusively by electrostatic forces, and is dependent on the ionization states of both lipid A and the drug, Under conditions when lipid A is highly disaggregated, pentamidine binds specifically to bis-phosphoryl- but not to monophosphoryl-lipid A indicating that both phosphate groups of lipid A are necessary for electrostatic interactions by the terminal amidininium groups of the drug, Based on these data, a structural model is proposed for the pentamidine-lipid A complex, which may be of value in designing endotoxin antagonists from first principles.

Structural studies of analgesics and their interactions. XII. Structure and interactions of anti-inflammatory fenamates. A concerted crystallographic and theoretical conformational study

A theoretical conformational analysis of fenamates, which are N-arylated derivatives of anthranilic acid or 2-aminonicotinic acid with different substituents on the aryl (phenyl) group, is reported. The analysis of these analgesics, which are believed to act through the inhibition of prostaglandin biosynthesis, was carried out using semi-empirical potential functions. The results and available crystallographic observations have been critically examined in terms of their relevance to drug action. Crystallographic studies of these drugs and their complexes have revealed that the fenamate molecules share a striking invariant feature, namely, the sixmembered ring bearing the carboxyl group is coplanar with the carboxyl group and the bridging imino group,the coplanarity being stabilized by resonance interactions and an internal hydrogen bond between the imino and carboxyl groups. The results of the theoretical analysis provide a conformational rationale for the observed invariant coplanarity. The second sixmembered ring, which provides hydrophobicity in a substantial part of the molecule, has limited conformational flexibility in meclofenamic, mefenamic and flufenamic acids. Comparison of the conformational energy maps of these acids shows that they could all assume the same conformation when bound to the relevant enzyme. The present study provides a structural explanation for the difference in the activity of niflumic acid, which can assume a conformation in which the whole molecule is nearly planar. The main role of the carboxyl group appears to be to provide a site for intermolecular interactions in addition to helping in stabilizing the invariant coplanar feature and providing hydrophilicity at one end of the molecule. The fenamates thus provide a good example of conformation- dependent molecular asymmetry.

Interactions of Antithyroid Drugs and Their Analogues with Halogens and their Biological Implications

The selenium analogue of antithyroid drug methimazole (MSeI) reacts with molecular bromine to produce two different types of novel complexes depending upon the molar ratio of MSeI to Br-2 in the reaction medium: Dicationic diselenide complex with two Br- ions as counterions is produced in the reaction of MSeI with 0.5 equiv of Br-2 (MSeI/Br-2, 1.0:0.5), whereas a stable 10-Se-3 hypervalent ``T-shaped'' complex featuring a linear Br-Se-Br moiety was produced when MSeI was treated with Br-2 in an equimolar ratio (MSeI/Br-2, 1.0:1.0). A substitution at the free N-H group in MSeI alters its reactivity toward iodine/bromine. For example, the N,N-disubstituted selones exclusively produce the corresponding 10-Se-3 hypervalent ``T-shaped'' complexes in the reaction with I-2. In the presence of the lectoperoxidase/H2O2/I- system, N,N-dimethylimidazole-2-selone produces the corresponding dicationic diselenide with two I- counterions as the final metabolite. The formation of ionic species in these reactions is confirmed by single crystal X-ray diffraction studies and in some cases by Fourier transform-Raman spectroscopic investigations.

Weak interactions involving organic fluorine: analysis of structural motifs in Flunazirine and Haloperidol

The crystal structure of Flunazirine, an anticonvulsant drug, is analyzed in terms of intermolecular interactions involving fluorine. The structure displays motifs formed by only weak interactions C–H⋯F and C–H⋯π. The motifs thus generated show cavities, which could serve as hosts for complexation. The structure of Flunazirine displays cavities formed by C–H⋯F and C–H⋯π interactions. Haloperidol, an antipsychotic drug, shows F⋯F interactions in the crystalline lattice in lieu of Cl⋯Cl interactions. However, strong O–H⋯N interactions dominate packing. The salient features of the two structures in terms of intermolecular interactions reveal, even though organic fluorine has lower tendency to engage in hydrogen bonding and F⋯F interactions, these interactions could play a significant role in the design of molecular assemblies via crystal engineering.

Crowd Sourcing a New Paradigm for Interactome Driven Drug Target Identification in Mycobacterium tuberculosis

A decade since the availability of Mycobacterium tuberculosis (Mtb) genome sequence, no promising drug has seen the light of the day. This not only indicates the challenges in discovering new drugs but also suggests a gap in our current understanding of Mtb biology. We attempt to bridge this gap by carrying out extensive re-annotation and constructing a systems level protein interaction map of Mtb with an objective of finding novel drug target candidates. Towards this, we synergized crowd sourcing and social networking methods through an initiative `Connect to Decode' (C2D) to generate the first and largest manually curated interactome of Mtb termed `3interactome pathway' (IPW), encompassing a total of 1434 proteins connected through 2575 functional relationships. Interactions leading to gene regulation, signal transduction, metabolism, structural complex formation have been catalogued. In the process, we have functionally annotated 87% of the Mtb genome in context of gene products. We further combine IPW with STRING based network to report central proteins, which may be assessed as potential drug targets for development of drugs with least possible side effects. The fact that five of the 17 predicted drug targets are already experimentally validated either genetically or biochemically lends credence to our unique approach.

Rationalization and prediction of drug resistant mutations in targets for clinical anti-tubercular drugs

Resistance to therapy limits the effectiveness of drug treatment in many diseases. Drug resistance can be considered as a successful outcome of the bacterial struggle to survive in the hostile environment of a drug-exposed cell. An important mechanism by which bacteria acquire drug resistance is through mutations in the drug target. Drug resistant strains (multi-drug resistant and extensively drug resistant) of Mycobacterium tuberculosis are being identified at alarming rates, increasing the global burden of tuberculosis. An understanding of the nature of mutations in different drug targets and how they achieve resistance is therefore important. An objective of this study is to first decipher sequence as well as structural bases for the observed resistance in known drug resistant mutants and then to predict positions in each target that are more prone to acquiring drug resistant mutations. A curated database containing hundreds of mutations in the 38 drug targets of nine major clinical drugs, associated with resistance is studied here. Mutations have been classified into those that occur in the binding site itself, those that occur in residues interacting with the binding site and those that occur in outer zones. Structural models of the wild type and mutant forms of the target proteins have been analysed to seek explanations for reduction in drug binding. Stability analysis of an entire array of 19 mutations at each of the residues for each target has been computed using structural models. Conservation indices of individual residues, binding sites and whole proteins are computed based on sequence conservation analysis of the target proteins. The analyses lead to insights about which positions in the polypeptide chain have a higher propensity to acquire drug resistant mutations. Thus critical insights can be obtained about the effect of mutations on drug binding, in terms of which amino acid positions and therefore which interactions should not be heavily relied upon, which in turn can be translated into guidelines for modifying the existing drugs as well as for designing new drugs. The methodology can serve as a general framework to study drug resistant mutants in other micro-organisms as well.

Network approaches to drug discovery

Introduction: Advances in genomics technologies are providing a very large amount of data on genome-wide gene expression profiles, protein molecules and their interactions with other macromolecules and metabolites. Molecular interaction networks provide a useful way to capture this complex data and comprehend it. Networks are beginning to be used in drug discovery, in many steps of the modern discovery pipeline, with large-scale molecular networks being particularly useful for the understanding of the molecular basis of the disease. Areas covered: The authors discuss network approaches used for drug target discovery and lead identification in the drug discovery pipeline. By reconstructing networks of targets, drugs and drug candidates as well as gene expression profiles under normal and disease conditions, the paper illustrates how it is possible to find relationships between different diseases, find biomarkers, explore drug repurposing and study emergence of drug resistance. Furthermore, the authors also look at networks which address particular important aspects such as off-target effects, combination-targets, mechanism of drug action and drug safety. Expert opinion: The network approach represents another paradigm shift in drug discovery science. A network approach provides a fresh perspective of understanding important proteins in the context of their cellular environments, providing a rational basis for deriving useful strategies in drug design. Besides drug target identification and inferring mechanism of action, networks enable us to address new ideas that could prove to be extremely useful for new drug discovery, such as drug repositioning, drug synergy, polypharmacology and personalized medicine.